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Galactic weakly interacting massive particles (WIMPs) may scatter off solar nuclei to orbits 
gravitationally bound to the Sun. Once bound, the WIMPs continue to lose energy by repeated 
scatters in the Sun, eventually leading to complete entrapment in the solar interior. While the 
density of the bound population is highest at the center of the Sun, the only observable signature 
of WIMP annihilations inside the Sun is neutrinos. It has been previously suggested that although 
the density of WIMPs just outside the Sun is lower than deep inside, gamma rays from WIMP 
annihilation just outside the surface of the Sun, in the so called WIMP halo around the Sun, may 
be more easily detected. We here revisit this problem using detailed Monte Carlo simulations and 
detailed composition and structure information about the Sun to estimate the size of the gamma- 
ray flux. Compared to earlier simpler estimates, we find that the gamma-ray flux from WIMP 
annihilations in the solar WIMP halo would be negligible; no current or planned detectors would 
be able to detect this flux. 



I. INTRODUCTION 

The observational evidence for dark matter is overwhelming (see e.g. [l|), but in spite of that, direct evidence for 
dark matter particles has not yet been found. One of the most popular classes of dark matter models is the Weakly 
Interacting Massive Particle (WIMP) model. WIMPs are heavy particles with weak interaction strengths, and are 
popular because they are naturally thermally produced at right abundances in the early Universe (see e.g. (2|). There 
are many ways to search for WIMPs, e.g. via direct detection, or via indirect detections, like neutrinos from the Sun 
[3, 0] or the Earth , gamma-rays from annihilations in the Galactic halo or charged cosmic-rays from annihilations 
in the Galactic halo (for a review, see [1]). The neutrino signal from the Sun arises from WIMPs scattering in the Sun, 
eventually sinking to the core and annihilating, producing many standard model particles. Of these, only neutrinos 
escape the Sun and the neutrinos can be searched for with neutrino telescopes like e.g. IceCube Q. 

We will here focus on another idea, proposed by Strausz [gj. The idea is that while WIMPs are being captured 
by the Sun, they will move on bound orbits outside the Sun, making up a halo of WIMPs around it. Pair-wise 
annihilations of WIMPs in this halo could occur and (if the density would be high enough) , produce a detectable flux 
of e.g. gamma-rays. This type of gamma-ray flux has been searched for by the Milagro [lO| detector. 

Strausz came to the conclusion that this gamma-ray signal would be measurable or that if no such signal is 
seen, one would be able to constrain parameters in the WIMP model. However, Strausz' calculation imposes many 
approximations and in a later preprint by Hooper jlli] , it was found that the flux of gamma-rays is much smaller than 
Strausz' estimates and a detection of this gamma-ray signal would require unrealistically large telescope areas. In 
both these calculations many simplifying approximations are made; for example radial orbits are assumed, the WIMPs 
real orbits are not followed with realistic solar models. Another calculation was performed by Fleysher p^ . finding 
even higher gamma-ray fluxes than in Strausz' calculation. However, many simplifying assumptions were made in all 
of these calculations, making the results unreliable. 

To clarify the discrepancies between these earlier estimates, we have here performed a much more thorough calcu- 
lation of the WIMP density around the Sun via detailed Monte Carlo simulations. We let WIMPs from the Galactic 
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FIG. 1: An example of a WIMP being captured by the Sun. A WIMP from the Milky Way halo population scatters at point 
1 as it passes through the Sun. If enough energy is lost the WIMP ends up on a bound orbit which it travels in, typically for 
many solar passages, until it scatters again in point 2 to an orbit with lower energy. In this example the WIMP has, after the 
scatter in point 3, lost enough energy to become completely trapped inside the Sun and is no longer of any interest to us. Note 
that in the figure the orbits are entirely elliptical, whereas in reality (and in our simulations), they are slightly non-elliptical 
inside the Sun. 

halo scatter in the Sun and follow the WIMPs on their orbits (taking full account of the non-Keplerian nature of the 
orbits inside the Sun), letting them interact repeatedly with solar nuclei until they are fully trapped inside the Sun. 
We calculate how much each WIMP we simulate contributes to the WIMP density around the Sun. By summing 
up the contributions from the entire WIMP population, we can estimate the total WIMP density around the Sun. 
We finally calculate the gamma-ray signal from WIMP annihilations in this halo, which can be compared with the 
Milagro searches for a gamma-ray signal of this kind [IQ] . We will also calculate predicted fluxes for future gamma-ray 
detectors and compare with expected backgrounds. 

II. THE WIMP CAPTURE PROCESS 

As Milky Way WIMPs pass through the Sun, some will scatter off solar nuclei and some of those will lose enough 
energy to become gravitationally bound to the Sun's potential well. The bound WIMPs will orbit the centre of the 
Sun and eventually scatter again in the Sun, losing energy in each scatter. Eventually enough energy is lost for 
the WIMP to end up on an orbit which lies completely inside the Sun (at which point we will called it completely 
trapped). The WIMP capture process is illustrated in Fig.[T] 

Due to the smallness of the scatter cross section, a bound WIMP typically survives many passages in the Sun before 
it scatters again. Also, as the WIMP is typically much heavier than the nuclei it scatters off the energy loss in each 
scatter tend to be quite small. Because of this the WIMPs often have to scatter many times before enough energy is 
lost for the WIMP to become completely trapped inside the Sun. In this way a captured WIMP can contribute to the 
WIMP density around the Sun for a long time even though a small energy loss in the first scatter can imply that the 
first orbit stretches very far out from the Sun. The effects of small energy losses in each scatter is most pronounced 
if hydrogen is an important target for WIMP scattering in the Sun. 

To get good accuracy of our predictions, we follow the WIMPs through the Sun on their actual orbits, taking the 
actual solar environment (as a function of radius), and details of the scattering process into account. 

All these intermediate orbits contribute to the WIMP density around the Sun, building up a WIMP halo around 
it (we will refer to this population of WIMPs as the solar WIMP halo, or just WIMP halo, not to be confused with 
the galactic WIMP halo, or galactic halo for short). WIMP annihilations within this halo give rise to hig h-energy 
gamma-rays, where the dominant background is cosmic-ray interactions with the solar chromosphere [l3[ and up- 
scattering of solar photons via inverse- Compton scattering from cosmic-ray electrons p^ - flGj . Also, the Sun is opaque 
to gamma-rays and hence shields against the diffuse gamma-ray background. 
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III. BUILD-UP OF THE SOLAR WIMP HALO 
A. The WIMPs' first scatter in tiie Sun 



Following the notation of Gould the total number of WIMPs that scatter in a shell of a given radius in the 
Sun per unit time and velocity is given by 

47rr w , (1) 



drdu 



where f{u) is the velocity distribution of the WIMPs in the Milky Way halo, in the Sun's reference system and 
infinitely far away, m is the mass of the nucleus in the Sun which the WIMP scatters off, a is the scatter cross section 
on that nucleus, u is the WIMP velocity at infinity, v is the escape velocity at the shell, w is the velocity of the WIMP 
in the shell, r is the radius of the shell and p(r) is the density of the element with nucleus mass m at radius r. We 
will later denote the WIMP mass with M. Eq. ([1]) should be summed over all elements in the Sun. We use the solar 
model of [18] with heavy element abundances from [l^ as implemented in DarkSUSY [2C|]. In total we include the 16 
most important elements in the Sun in our calculation. 

Assuming the velocities of the WIMPs in the Milky Way halo to be Maxwell-Boltzmann distributed, the velocity 
distribution in the Sun's reference frame is [l7| 

f{u) _ [6 1 / 3w|\ f 3u^\ . f3uvQ\ 



U \ TT VqV 



du = nwW r exp --rrr exp -— sinh du, (2) 



where Vq is the Sun's velocity orbiting the Milky Way centre and v is the WIMP three-dimensional velocity dispersion 
in the Galactic halo, here taken to be w = 270 km/s. nw — Pw/M is the number density of WIMPs in our region of 
the Milky Way halo, where we take pvv = 0.3 GeV cm^'^. 

For further use the WIMPs' energy distribution will be more useful than the velocity distribution; Eq. ([2]) can be 
rewritten as 



£ 



[3 1 / 3w?,\ / 3£\ . , f3vQ r—\ 1 ,^ 
at = nw\ exp — -— r exp — ^ smh -r:T-v2t —= at, 

V TT VqV \ 2v^ J \ V J V « J V£ 



(3) 



where £ is the reduced energy, defined as the WIMP's energy (kinetic plus potential) divided by the WIMP mass (for 
a captured WIMP £ < 0). Analogously is the reduced kinetic energy and £p is the reduced potential energy at 
some radius. 

Rewriting Eq. ([T]) now yields the number of WIMPs which scatter in a shell in the Sun per unit time and reduced 



energy ^ 



d£ dr m £ 

expressed using the reduced energy the WIMP has before it scatters 



1. The WIMPs' energy distribuUon after their first scatter 

When a WIMP interacts with a nucleus it scatters from the reduced kinetic energy before the scatter, fk, to the 
reduced kinetic energy after the scatter, fk, which, as for general elastic scatter, fulfils the relation 



M + m 



fk < < £k (5) 



1 Note that 

liin M = but lim ^ 

ii->0 U £-»0 S 

and hence the number of scatters per unit time and velocity vanishes as m — > while the number of scatters per unit time end energy 
does not vanish as £ — > 0. 
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with all £k in this interval being equally probable. Note that it is the energy loss and not the loss in velocity that is 
uniformly distributed. 

This means that a WIMP with the reduced energy £y can only have scattered from reduced energies fulfilling 



< < fk, (6) 

\M — m J 

which written in terms of total reduced energy E — E]^ + E-p (and E — E]^+ E-p, with E-p — Ep) becomes 

E<E<[ E - — -Ep. 7 

\M ~ m J [M — my 

Note that the reduced energies are greater than zero for particles not bound to the solar system, hence £ > for 
the WIMPs in the Milky Way halo, whereas E can be both negative or positive, depending on the energy lost in the 
scatter. 

The number of WIMPs which scatter to a certain reduced energy, in a given shell in the Sun, per unit time is 
then 



d^r _ ^^^^f^jm^Mf r\imei^,^,E. 

m Mm E ^ ' 



dE dr 



(8) 



[M + mV- AMm ^ 
where u.L. — — E — -r— r^c, 



\M-mJ (M-m) 



which is the integral of ^ with integration limits given by (O. Here 9 is the Heaviside function, which is required 
since E can be negative while E is always positive. 

The above expression has been normalised by multiplying by the reciprocal of the interval length in ([5]) 

1 (M + m)2 



fk 4Mm 



2. Form factor suppression of the scattering cross section 

For scatter off nuclei larger than the proton the cross section depends on the energy loss in the scatter. For scatters 
with low momentum transfer the wave functions of the different nucleons in the nuclei may be viewed as coherent, 
making the cross section increase as the square of the number of nucleons in the nucleus, A^. For higher momentum 
transfer this coherence breaks down. This suppression is taken into account by introducing an exponential form factor, 
suppressing high-energy transfers, i.e. high AE [l^ 

expf-^) (10) 



with 



£a = ^ and R • 



0.91 (^) +0.3 



X 10~^^ m, (11) 



where R is an estimate of the radius of the nucleus [21 1 



Introducing a form factor, Eq. can be written as 



ST _ ,„,£W r' M) e.p eiew. (12) 



dEdr m Mm Ji E \ Eq 

fM + mV~ AMm ^ 
where u.l. = — E - —- -rrEp 

where a refers to the total low-energy cross section (i.e. without the form factor suppression). Since the Sun contains 
many different elements, p(r) in (jl2[) should be viewed as the density of a specific element and then summed over all 
relevant elements in the Sun, with the form factor removed for hydrogen. 
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FIG. 2: Plot of Eq. (|12p for the case where spin-dependent scatter, ctsd = 10~^ pb, and spin-independent scatter, ctsi ~ 10~^ 
pb, respectively being dominant. The graphs show the distribution of scattered WIMPs after their first scatter in the Sun, 
expressed as the number of scattered WIMPs per unit time, unit reduced energy after the scatter and unit radius in the Sun. 
The different lines correspond to different radii in the Sun where the WIMPs scatter, showing the shells with radius 10%, 20%, 
. . . , 60% of the Sun's radius, the graphs hence illustrates surfaces. The smaller the radius of the illustrated shell, the further 
to the left the graph has non-zero values, i.e. the smaller the possible reduced energy is after the scatter. Also indicated are 
the reduced energies corresponding to orbits completely trapped in the interior of the Sun. The WIMP mass is here set to 100 
GeV. 



Regarding the cross sections, WIMPs can typically scatter via spin-independent and spin-dependent scatterings. 
The spin-independent scatterings occur on all elements in the Sun, whereas the spin-dependent only occur for nuclei 
with spin. The only element with spin with non- negligible abundance in the Sun is hydrogen. In Eq. ()12|) . the cross 
section a refers to the total cross section on a given element in the Sun (without the form- factor suppression). Usually 
one relates this to the cross section on a single proton, which we will denote csd and asi for the spin-dependent 
and spin- independent cases respectively. For the spin-dependent case a = usd (as only scatterings on hydrogen are 
relevant), whereas for the spin- independent case, the cross sections are related via 

^ \AI + mJ \M + mpJ ^ ' 

where rup is the mass of the proton and A is the atomic number of the target nucleus. 

Eq. ([T^ describes the distribution of WIMPs after their first scatter in the Sun, illustrated in Fig. [51 Hydrogen 
is the only significant element in the Sun which couples axially to the WIMP so if spin-dependent scatter, (Jsd, 
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dominates the WIMPs will scatter effectively only off hydrogen. If spin-independent interactions dominate, scatter 
off hydrogen is insignificant and the WIMPs will essentially only scatter off helium and heavier elements in the Sun. 
Heavier elements can play an important role even though their abundance is small since the scatter cross section 
increases for large nuclei and because of the increase of phase space accessible in the scatter. As the WIMPs scatter 
off nuclei of different masses in the two scenarios the dynamics is quite different. Again, note that both crgD and crgi 
refer to the scattering cross sections off protons (as these are the quantities usually used to quantify the size of the 
scattering cross sections). We will use optimistic values of tisi = 10~^ pb and ctsd = 10""^ pb, but as we will see in 
section HV Al our results do not depend on the values chosen for the cross sections (within reasonable limits). 

In the upper graph in Fig. [2] the WIMPs scatter only off hydrogen which makes the energy loss in the scatter quite 
small, making the Sun only capable of capturing the low energy tail of the Milky Way WIMP population. This is 
seen in the graph in that the majority of the WIMPs scatter to an energy greater than zero. As we will discuss later, 
the WIMPs have to scatter many times before being fully trapped inside the Sun, allowing the bound WIMPs to give 
larger contributions to the WIMP density outside the Sun. 

In the lower graph in Fig. [2] the spin-independent scatter dominate making the WIMPs scatter mostly off heavier 
nuclei and hence, on average, lose a lot more energy in the scatter. This makes the WIMPs much more likely to be 
captured in the Sun's potential well but as seen in the graph some will immediately disappear inside the Sun since 
a WIMP needs a reduced energy of at least £ = ——^^ = —1.9 x 10^^ m^/s^ to be energetic enough to reach the 

region outside of the Sun. For circular orbits, £ = —3^^ = —0.95 x 10^^ m^/s^ is the limiting energy for the orbits 
being fully trapped inside the Sun. These two regions are shown in the lower panel in Fig. [2l As the orbit of a bound 
WIMP must contain the point of last scatter (up to the spherical symmetry) , the WIMP orbits after first scatter are 
necessarily radial rather than circular. With this in mind we see that it is not very likely that the WIMP loses enough 
energy to be completely trapped already in the first scatter. For heavier masses than the M = 100 GeV shown in 
Fig. [21 it is even harder to fully capture the WIMP with only one scatter. For spin-dependent scatters (upper panel), 
not enough energy can be lost for the WIMPs to be completely captured in the Sun with only one scatter. 

The pointiness of the graphs in Fig. [2] comes from the discontinuity of the integrand which is related to the footnote 
associated to Eq. (g]). The function ^^9{£) jumps from zero to its maximum value at £ — 0. 

WIMPs which scatter close to the centre of the Sun have higher velocities in the scatter and hence normally lose a 
larger fraction of their total energy in the scatter. This is the reason for inner shell curves in Fig. [2] being non-zero 
for lower values of £ . 

Integrating Eq. (fT2|) over r and £, for £ < 0, gives the total capture rate of WIMPs by the Sun. The total capture 
rate is of interest in computing the neutrino flux from annihilating WIMPs in the centre of the Sun. This can also be 
calculated using the DarkSUSY program [20] and is in excellent agreement with our results. 



3. The scatter angle after the WIMPs ' first scatter 



To determine the orbit of the captured WIMP one, besides the energy distribution and the place of scatter, also 
needs the scatter angle, a, defined as the angle between the path of the scattered WIMP and the radial direction 
from the centre of the Sun. 

In this work spherical symmetry is assumed; the Sun is spherically symmetric and isotropy in the incoming WIMPs 
is assumed. The small scatter probability for a WIMP passing through the Sun makes it equally likely for the WIMP 
to scatter on its way out as on its way in to the Sun. Because of this the scattered WIMPs inherit the initial WIMP 
isotropy and all the directions of the just scattered WIMPs are equally probable. If one is to be very careful, the Sun's 
velocity and gravitational infiuence could give rise to some small anisotropy. As the orbits are not purely elliptical 
inside the Sun, they will precess, which will somewhat smoothen out anisotropics. Gravitational interactions with the 
planets will also have a similar effect on orbits reaching further out from the Sun. 

All scatter directions being equally probable after the first scatter gives the probability distribution for a 

27rsin(a)da 1 . , x, , „ ,^ ,^ 
— - smlalda, where < a < tt. 14 



B. The orbits and lifetimes of the captured WIMPs 



The reduced energy, £, and reduced angular momentum, J7 (defined as the angular momentum divided by the 
WIMP mass M), of a captured WIMP are given by 

£ = y+^pM, (15) 
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J^rv\suia\ = ry^2(£-£p(r))|sina|, (16) 

where r is the radius where the scatter took place and v the velocity the WIMP scattered to. The reduced potential 
energy at the scatter, Ep, is given by a numerical function for the gravitational potential inside the Sun. 

The energy and angular momentum together fully specify the orbit, up to the spherical symmetry. It will be useful 
for us to find the minimal and maximal distances of the orbit to the centre of the Sun, r^m^ and r,„ax- These two 
points are the only points in the orbit fulfilling J — r^2(E — E^). We are here only interested in WIMPs which 
spend time both inside and outside the Sun, i.e. WIMPs on orbits fulfilling rmin < -R© < '"max- As the WIMPs were 
captured by scatters in the Sun their future orbits will always be partly inside the Sun unless disturbed in some way, 
e.g. by the planets; this is not taken into account in these calculations. Outside the Sun the gravitational potential is 
simply 5p = —GMqIt^ allowing us to simply solve for r„iax 

rmax = -^ + y(^^j +2£- 

Inside the Sun the gravitational potential is not so simple and hence rmin can not be found in such a simple way as 
^max- The orbit everywhere fulfils JT^^ ^ {vv)^ — 1t^{E 6'p(r)) with equality only for v equal to Tmin or rmax- 

To 

find rniin, we can use that > 2r^{E — Ep{r)) for r < rmin and perform a numerical binary search for rmin. 



1. The density contribution to the WIMP halo before the WIMP scatters again. 

To find the density of WIMPs around the Sun one needs to know what fraction of its time an orbiting WIMP 
spends at different distances from the Sun. The density contribution from a WIMP at a certain radial interval is 
proportional to 2St/T, where St is the time it takes the WIMP to traverse the shell and T is the period of the orbit. 
The factor of two is added as in completing an orbit the WIMP passes all the radii in the orbit twice. 

We will later look closer at the radial velocities, and from Eq. (|23l) one finds, for the Keplerian potential outside 
the Sun, that the time it takes to traverse a shell of thickness |dr| is 

dt^ , =|dr|. (18) 

V2(GMo + Er}r - 

This expression diverges for r — rmin and r — rmax but the integral converges everywhere. Since the velocity is lower 
in the outer regions of the orbit the WIMP also spends most of its time there. 

To reduce noise in extracting the density from the Monte Carlo the density contribution from a bound WIMP is 
given as the integral of the density contribution over the (small) radii interval which is the bin width in the density 
plots. This is easily done as the integral of Eq. (|18|) can be found in closed form; using Mathematica a primitive 
function of the right hand side of Eq. is given by 

- . arctan ( '2Er + GM^ ^ -^2r{GMr^ + Er) - J\ (19) 



(-2f)3/2 \^r2E^2r{GMQ^Er)-J' 



^ + iV2r(GMo+£:r)-^2_ 



As discussed above a factor of 2, compared to Eq. ([T8|. has also been added to the primitive function since an orbiting 
WIMP passes a radial interval twice or not at all. 

For the orbits investigated here the WIMP typically spends most of its time outside the Sun. The total time, T, it 
takes for the WIMP to complete a full revolution in an elliptical orbit is 

^^GM^ (20) 

The density contribution also depends on how many passages the WIMP survives before it scatters again. For each 
solar passage there is some probability for the WIMP to scatter in the Sun, reducing the average density contribution 
from a given bound orbit for each solar passage. The probability for the WIMP to scatter twice in one solar passage 
is immensely small and can be neglected. 

The probability for the bound WIMP to survive n solar passages without scattering in the Sun is = (1 — P)" 
with P being the probability for the WIMP to scatter in one passage. The total density contribution by a WIMP in 
a given orbit, integrated over time, is, on average 

9 „x l-g"+^ l-(l-P)"+i , , 

1 — g P 
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where po is the density contribution from the WIMPs first revolution in the given orbit. The finite age, of approx- 
imately 4.5 billion years, of the Sun gives a cut in the sum since no orbits can be older than this. The maximum 
number of revolutions a WIMP can have fulfilled in a given orbit is n = (4.5 x 10^ yr)/T, with T being the orbit time, 
as in expression (^0)) . 

In practice, the finite age of the Sun is important only for very low scatter cross sections. Then also the orbits 
which these long lived WIMPs scatter to should be removed. This was taken into account in our runs with very low 
scatter cross section. 



2. The scatter probability for a bound WIMP passing through the Sun 



The probability for collision when a WIMP passes a shell in the Sun is given by 

p(r) „ dr 

m \Vr\ 



(22) 



where p and m refer to the density and mass, respectively, of the target nucleus. The above expression is then summed 
over the relevant elements. Also, the scatter cross section can depend on the energy loss in the scatter, which depends 
on the WIMP velocity and hence the radius of the scatter. The above expression comes from the number of scatters 
per unit time {anv) times the time the WIMP spends in a given shell taken twice, since the WIMP intersects a shell 
twice or not at all in one solar passage. 

From V and for a particle one easily determines the radial velocity component 



V 

R 



Hence, the probability for the bound WIMP to scatter in the Sun in one solar passage is 



Rp, 



P = - / ap{r) 



: dr, 



(23) 



(24) 



where ^ 2£ — 2£p(r) and r„iin is easily determined numerically, as previously discussed. If the WIMP scatters off 



hydrogen, the cross section does not suffer a form- factor suppression and it is straightforward to evaluate Eq. ([M]) . 

Accounting for the form factor suppression of the cross section needs some further treatment. For a bound WIMP 
passing the Sun once in an orbit with reduced energy, £, and reduced angular momentum, J, the probability to 
scatter off a given element i in a (possible) shell in the Sun, given by the radius r, to a certain, possible reduced 
energy £ is given by 



drd£ 



1 (M 



~ P^ir) 2 



(25) 



£k 4Mmi 

where the rightmost part is simply Eq. (|22l) . The first part of Eq. (P5|) is the factor given in Eq. ^ times the form 
factor suppression. 

The probability for the orbiting WIMP to scatter to any possible energy in a given shell during one solar passage 
is given by the integral of Eq. (j25|) over the possible values of £ and then summed over the different elements in the 
Sun 



dP 
dr 



V 1 



2(o-sD + o-si) 



mi 



16 



^(M + m^., p,(r) 



2Mm, 



1 — exp 



4Afmj 
" [M + m,Y £^ 



(26) 



where v/vy is given by Eq. (1231) . The mass and density of element i are given by m^ and Pi{r) respectively, with 
i = 1 being hydrogen. The leftmost term refers to scattering off hydrogen, which is the only element here with 
spin-dependent scatter cross section to the WIMP (the spin-independent scattering off hydrogen is added to this term 
as well for clarity). (Tgj here refers to the total low-energy spin- independent cross section on element i, and is given 
by a in Eq. P^ . 

The probability for the captured, bound WIMP to scatter during one solar passage in its orbit is then given by the 
integral of the above expression over the radii in the Sun that the WIMP passes through in its orbit. 
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FIG. 3: The line density n around the Sun of WIMPs that have scattered exactly once in the Sun. The black curve is for 
spin-dependent scatters dominating, using asD = 10~^ pb, asi = and is generated by a calculation of 1 million scattered 
WIMPs. The red curve is for spin-independent scatters dominating, using asi = 10~^ pb, aso = and is the result of 100 000 
simulated scattered WIMPs. For both curves the WIMP mass is M = 100 GeV. The Earth is approximately 215 solar radii 
away from the Sun. 



3. The density contribution of the WIMPs which have only scattered once 

The analysis is now sufficient to calculate the density around the Sun of WIMPs which have scattered exactly once 
in the Sun. This is done through making a Monte Carlo simulation, simulating a large number of captured WIMPs 
according to the energy distribution after their first scatter, given by Eq. p^ . This is described in the first part of 
section UlI Dl which describes the structure of the final full Monte Carlo. The resulting density contribution of the 
WIMPs which have scattered only once is shown in Fig. [3l This graph shows the line number density, n, of WIMPs 
around the Sun. The line number density is the number of WIMPs in a spherical shell of a given radius r per radial 
interval and it is related to the three-dimensional spatial density N (number of WIMPs per volume element) as 

- 3 (27) 

The line number density in Fig. [3] is the contribution to the line number density from the WIMPs that have scattered 
once, but not yet scattered a second time, and become gravitationally bound to the Sun. As they scatter more times, 
they will move in closer to the Sun and the contribution will increase at lower radii. The reason the spin-dependent 
curve peaks at higher radii is kinematical as the WIMPs in this case have only scattered on hydrogen. The spin- 
independent curve includes scatterings on heavier elements, which causes the WIMPs to lose more energy in the 
scattering. 

The unevenness in Fig. [3] is due to the stochastic nature of Monte Carlo simulations. We have simulated a finite 
number of WIMP orbits to calculate the contribution to the line number density and as some simulated orbits can 
be quite stable (meaning that they have a long life time and hence a large contribution to the line number density), 
we can get large contributions from a few orbits. The peaks then occur at the turn-around radius rmax, where the 
WIMPs spend most of their time and these are the peaks we see in the figure. The more WIMP orbits we simulate, 
the more even the curves get, but given the fact that we will be most interested in the density close to the Sun (where 
the density is highest, and the curves are smoother), we do not have to worry about peaks at large radii. Also, the 
spin-dependent curve is smooth at small radii as in this case the WIMPs close to the Sun are never close to the 
turn-around radius Tmax after the first scattering. 
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Taking into account the fact that the orbits are not perfect elhpses inside the Sun is found to increase the calculated 
line number density by roughly 20%. This is because the true orbits do not come as far into the dense central regions 
in the Sun as purely elliptical orbits would, which increases the life time of the true orbits since the scatter probability 
decreases. 



C. Subsequent scatters in the Sun 

In general a WIMP scatters many times before its orbit becomes completely trapped inside the Sun. These 
intermediate orbits contribute to the WIMP density, especially in the region close to the Sun. 

The radius of the point of scatter for an already bound WIMP and the type of nucleus involved can in the Monte 
Carlo be extracted from Eq. (f26| . since this encodes the distribution of the scatter parameters for a given orbit. 
Knowing the place of scatter and the elements involved, the distribution of the energy loss in the scatter is found by 
expression pS]). 

To find the WIMP's angular momentum after the scatter, J, requires some further treatment. Once the properties 
of the new orbit are determined the contribution to the halo density can be calculated in the same way as for the 
WIMPs which have only scattered once. 

Using conservation of energy and momentum in the scatter and the general relation — p ■ p ~ 2A4Ey^ one can 
derive the relation 

p-p = (M ~m)Ei, + {M + m)Ei„ (28) 

where p and p are the momentum before and after the scatter, respectively. Dividing p and r in components parallel 
and r||) and orhogonal (p_L and rj^) to p, using ([SS)) and p • p = pp\\ gives 

^ (Af - m)E^ + (M + m)E^ ^2 ^2 

P|| = and p^^p (29) 

From the angular momentum before the scatter, J = |J| = r±p, we can derive r± and r|| since both p and J are 
known. 

The modulus of the angular momentum after the scatter can be evaluated as 

2 

= r'^p\ sin^ (j) + {r\\p±^ cos4>- r±p\\)'^. (30) 

The azimuthal angle (f> is arbitrary in any given scatter and is hence generated as a random number in the Monte 
Carlo simulation. 



D. Monte Carlo 

The actual calculation of the density of the Sun's WIMP halo is implemented through constructing a Monte Carlo 
simulating the capture process of a large number of WIMPs. The Monte Carlo is constructed as follows. 

1. The Monte Carlo starts from the distribution of WIMPs which have scattered once to bound orbits. This 
distribution of WIMP interaction points, energy after the scatter and scatter element (from the solar model) is 
given by Eq. this distribution is also shown in Fig.Hl Having picked a WIMP according to this distribution, 
the scatter angle is picked from Eq. (|14p , giving the final piece of information required to calculate the WIMPs 
angular momentum. 

2. The density contribution by the WIMP in the orbit, before it scatters again, is then calculated. The line density 
contribution in a radial interval for one revolution in the orbit is given by Eq. (jl9p . This is then multiplied with 
the average number of solar passages that a WIMP in such an orbit survives, which is given by Eq. (PT|) and 

Eq. dm). 

3. The shell in the Sun in which the WIMP scatters next time is picked, using the expression for the scatter 
probability in the passage of a shell, given by Eq. The type of element the WIMP scatters off in the shell 
is then chosen according to the elemental abundance in the Sun at that radius. 







f pj_ COS0\ 






p±sin<f) 
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4. The WIMPs energy loss in the scatter and its new angular momentum is determined, as described in section 

5. If the WIMP's new orbit is still energetic enough to be partly outside the Sun (rmax > Rq in Eq. ^T7\ ). the 
simulation continues with step [H continuing to add to the WIMP density around the Sun. 

6. If the WIMP has lost enough energy to be completely trapped inside the Sun it is no longer relevant here. The 
simulation of a new scattered WIMP starts from step [1] 

7. When the desired number of captured WIMPs are simulated the result is weighted with the capture rate, which 
is found by integrating Eq. over r and £ ioT £ < 0, to get the real density. 



IV. THE WIMP NUMBER DENSITY IN THE HALO 



The constructed Monte Carlo has, as described in the above section, simulated the capture process, giving the 
density of WIMPs around the Sun. The simulations have been made for different WIMP masses and for both spin- 
dependent and spin-independent cross section being dominant. In total we have simulated about 10^ WIMPs from 
first scatter to complete solar entrapment. The results of these simulations are shown in Fig. |4l The line number 
density can be well approximated by 



-0.48 

^SD ^ ^025.21-1.015. I ^-1^ (3^) 



2 / r 



-0.40 



^i)J _ -^QZd.tl-V.ZiSZx-U.Wbbx j m^ (32) 

with X = logio • ^^^^ 

The three-dimensional spatial WIMP number density, N, is then given from the line number density by Eq. (j27p . The 
fits given above are optimized for small radii, as these are the most important for calculating the flux from annihilation 
in the solar WIMP halo. Note that the densities we find here are only marginally larger than the background density 
of WIMPs in the Galactic halo (about a factor of ten larger at the solar surface), so the density enhancements are 
not very dramatic. 



A. The halo number density dependence on the cross section and WIMP mass 



The magnitude of the cross section, except for very low cross sections, is insignificant for the WIMP density around 
the Sun. This is due to the capture rate of Galactic WIMPs in the Sun being proportional to the cross section, 
while the lifetime of a WIMP orbit is proportional to the inverse of the cross section. This cancels the cross section 
dependence if equilibrium has been reached, which happens if the Sun's capture rate of WIMPs from the Galactic 
halo equal the rate at which captured WIMPs become fully trapped inside the Sun. For very low scatter cross sections 
this process becomes so slow for a substantial fraction of the WIMPs that the age of the Sun is not long enough for 
equilibrium to have been reached, reducing the WIMP density around the Sun for very low scatter cross sections. 
The internal relation between the spin-dependent and spin-independent cross section is however important, since this 
determines whether or not the WIMPs scatter mostly off hydrogen or mostly off heavier elements. 

The WIMP mass density {p = nM) around the Sun is seen to not be all that sensitive to the WIMP mass; Fig. 2] 
shows the WIMP (line) number density which depends on WIMP mass, as does the WIMP number density in the 
Milky Way WIMP halo. For heavy WIMPs the mass difference between the WIMP and the target nucleus is larger, 
making it more unlikely for a Galactic halo WIMP to lose enough energy in a scatter to become gravitationally bound 
to the Sun. On the other hand a heavy WIMP also requires more scatters to lose enough energy to be completely 
trapped inside the Sun, allowing each captured WIMP to contribute to the dark matter density around the Sun for a 
longer time. These two effects somewhat cancel, making the WIMP mass density largely insensitive to WIMP mass. 
Note that as heavier WIMPs lose their energy more slowly, the WIMP density in this case becomes slightly more 
sensitive to very low scatter cross sections. 
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FIG. 4: a) The line number density of WIMPs around the Sun with spin-dependent cross section being dominant, i.e. WIMPs 
do only scatter off hydrogen. In this figure ctsd = 10~^ pb and asi = pb. The dashed curves are the fits described by 
expression (1311) . b) The line number density of WIMPs around the Sun with spin-independent cross section being dominant, 
i.e. WIMPs scatter off all elements in the Sun. In this figure ctsd = pb and ctsi = 10~^ pb. The dashed curves are the fits 
described by expression H32|l . 



V. GAMMA-RAY FLUXES AT EARTH 



A. Signal fluxes from WIMP annihilations in the solar halo 



The number of WIMP annihilations, Fa, for a given volume is in general given by 

Fa = J n'^dV. (34) 

For simplicity, let us assume a standard annihilation cross section of {<tv) = 3 x 10~^® cm'^s"^, that gives about the 
correct relic density of WIMPs. 

The flux of gamma radiation at Earth from a small source is 

= TA^iV^, (35) 

where / is the distance between the gamma-ray source and the telescope at Earth. The average number of gamma-rays 
created per WIMP annihilation, N^, requires further assumptions about the WIMP candidate investigated. We here 
assume that the WIMP annihilations create bb, which give reasonably high fluxes of gamma-rays. With DarkSUSY 
5.0.4 0, we have calculated that we then get about = 14.1,45.7 and 94.6 for M = 100, 1000 and 10000 GeV 
respectively. These are calculated above 1 GeV in gamma energy. We could also get monochromatic gamma-rays 
from annihilation to 77 and Zj or a harder spectrum from annihilation to e.g. T+r~, but (as we will see), given the 
low fluxes of gamma-rays we predict, it does not make any major difference to our conclusions including these other 
channels or not. 

Combining expressions (|34p and ([55)1 for our spherically symmetric WIMP distribution, one obtains the observed 
flux at Earth per solid angle 




which can be evaluated using = P -\- B? — 2RI cos{6), where R is the distance between the Earth and the Sun, 6 is 
the angle between the direction of observation and the direction towards the solar centre and r is the distance from 



13 




FIG. 5: The flux per solid angle from the Sun's WIMP halo as seen from Earth for M = 100 GeV, ctsd = 10"^ pb, osi = 
and A'^-, = 14.1. Plotted against the angle of sight 6, looking towards the centre of the Sun corresponds to 6 = 0. 



the Sun. Since gamma-rays can not pass through the Sun one gets a limit in I at the Sun's surface in the integration, 
when the line of sight intersects the Sun, giving 

l<Rcose~ ^JrI- sin^e when < arctan (^^^ • (37) 

The gamma-ray flux at Earth per solid angle, from Eq. is shown for a 100 GeV WIMP in Fig.[SJ When looking 
just outside the rim of the Sun the graph peaks since the distance probed by the telescope suddenly increases (as does 
the astrophysical background). The gamma-ray flux decreases fast with increasing 6 since the WIMP overdensity falls 
off fast with distance. Note that we have here not included the background density of WIMPs from the Galactic halo, 
that would actually dominate for the large distances away from the Sun. 

The total flux at Earth from these WIMP annihilations integrated out to some maximum angle, ^max, is 

<^>^^hav)N^ [ [ n^{r)smedlde. (38) 
4 Jo Jo 

In Table U we show this flux for some WIMP parameter configurations (integrated over all 9, as the end result is 
not very sensitive on the actual value of ^max)- One should also note that the WIMP annihilation rate in the solar 
halo is much lower than the capture rate and hence has no significant impact on the WIMP density in the solar halo. 

Given an experiment with an effective area A and live-time T, we can write the total number of events (total 
number of photons) as 

A^cvcnts = (39) 

where 



F = AT. (40) 

In reality though, the effective area A will depend (more or less strongly, depending on experiment) on energy, and a 
more accurate way of calculating F is to convolve the gamma-ray energy spectrum (where we use bb as our template 
spectrum) with the integrated (over time) effective areas, A'^g, 
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M (GeV) 


100 


100 


1000 


1000 


10 000 


10 000 


o-si (pb) 





10"^ 





10"^ 





lO"'' 


crSD (pb) 


10~^ 





10"^ 





10"^ 





A^T.(> 1 GeV) 


14.1 


14.1 


45.7 


45.7 


94.6 


94.6 




8.4 X 10"^^ 


1.8 X 10"^^ 


2.5 X 10"^° 


1.7 X 10"^" 


4.9 X 10"^^ 


4.0 X 10"^^ 


^Milagro (m^ s) 








2.2 X lO'' 


2.2 X 10^ 


3.6 X 10'' 


3.6 X 10^ 


A^Miiagro (events) 








5.6 X 10"^^ 


3.8 X 10"^^ 


1.8 X 10"" 


1.5 X 10"" 


-Fhawc (m^ s) 


2.1 X 10^ 


2.1 X lO'' 


1.6 X 10^ 


1.6 X lO'' 


7.9 X lO" 


7.9 X 10^° 


A^HAwc (events) 


1.7 X 10"^^ 


3.7 X 10"^^ 


3.9 X 10"" 


2.6 X 10"" 


4.0 X 10"" 


3.3 X 10"" 


-Fargo-ybj (m^ s) 


2.8 X 10^° 


2.8 X 10"' 


3.3 X 10" 


3.3 X 10" 


4.5 X 10" 


4.5 X 10" 


A''argo-ybj (events) 


2.4 X 10"** 


5.1 X 10"^ 


8.2 X 10"^ 


5.6 X 10"'' 


2.3 X lO"'-* 


1.9 X 10"^ 


Fpermi (m^ s) 


2.2 X 10^ 


2.2 X lO'^ 


2.3 X 10^ 


2.3 X lO'^ 


2.3 X 10'' 


2.3 X 10^ 


A^'pcrmi (events) 


1.9 X 10"" 


4.0 X 10"" 


5.7 X 10"^^ 


3.9 X 10"" 


1.1 X 10"" 


9.3 X lO"^'' 



TABLE I: A set of benchmark WIMP models and their corresponding gamma-ray fluxes above 1 GeV at Earth from WIMP 
annihilations in the solar halo. Also shown are the integrated effective areas convolved with the gamma-ray spectrum dN-,/dE 
for WIMP annihilations to bb, F. The total number of events expected in the Milagro search for gamma-rays from the Sun is 
also shown. For Milagro, the number of events are given for the specific search for gamma-rays from the Sun that has been 
carried out in For HAWC, ARGO-YBJ and Fermi, 5 years of observation of the Sun has been assumed with the effective 
areas taken from [2^ . [2^ and [Ij. The annihilation cross section is assumed to be {av) = 3 x 10"^'' cm""^ s~^. 

which is the expression we will use to calculate the number of events for some relevant gamma-ray experiments. 

Milagro has searched for gamma-rays from the Sun and with their estimated integrated effective areas A^^ for these 
searches [l(S|, based on 1165 hours of observation of the Sun, we can estimate the total number of events A^Miiagro these 
WIMP annihilations should have produced. In TableUwe give both F and the total number of events for Milagro. As 
can be seen, the total number of events are at best in the 10"^'^ range, an extremely small number of events, which 
of course would not have shown up in the Milagro searches at all. 

Aside from Milagro, there are other gamma-ray experiments {Fermi '25], HAWC ^26^, and ARGO-YBJ [23*1) that 
might be sensitive to a gamma-ray excess from the Sun. Air Cerenkov telescopes cannot look at the Sun as they can 
only observe when the sky is dark. Fermi has an effective area of slightly below 1 m^ and a field of view of about tt sr. 
We have integrated the effective area (integrated over field of view) for Fermi [23| , which for five years of operation 
gives a total number of events less than 10"^", far below any detection limit. HAWC on the other hand will be much 
larger than Milagro and Fermi and will have an effective area about 100 m^ at 100 GeV and about few times 10^ m^ 
at higher energies [l^ . We have calculated F for HAWC and the corresponding number of events in five years livetime 
and the number of events are at most in the 10"^^ range, i.e. very small. Even though HAWC is much bigger than 
Fermi, it has a rather high threshold of around 30 GeV and quite low effective areas for low energies, which makes 
the rates rather small for typical WIMP masses. ARGO-YBJ is designed to have a low energy threshold (sub-GeV) 
and the effective area increases rapidly with energy. At 1 TeV, the effective area is at most 10^ m^ [2^ (assuming 
the optimistic 1 multiplicity channel for scaler mode operation). We have calculated F also for ARGO-YBJ given 
the effective areas for this mode in [2^ . These and the corresponding number of events are also given in Table ID 
The number of events are at most in the 10"^ range for five years live-time. All these values for F and the expected 
number of events are given in Table ID The signals are really diminutive, even with the optimistic assumptions for 
current and future detectors used here. 



B. Backgrounds of gamma-ray fluxes from the Sun 

As the Sun does not shine at such high photon energies and also shields against the diffuse gamma-ray background, 
one might naively assume that this signal has essentially no background. However, secondary cosmic-ray processes 
produce gamma-rays from around the Sun. Cosmic-ray protons interacting with the Sun's chromosphere [13| . produce 
a gamma-ray flux of about 10"^ m~^ s~^ above 1 GeV. 

Solar photons can also be up-scattered in energy by cosmic-ray electrons p^ - [l6| . which produce an even larger 
flux of gamma-rays than the proton interactions. This gamma-ray flux has also been measured by EGRET data [27i |. 
where a flux of about 3 x 10"^ m~^ s~^ in the energy range 100-300 MeV has been seen. Also Fermi has observed 
this gamma flux in the energy range 100 MeV-10 GeV and the preliminary results agree well with the predicted 
background ^2§\. One should also note that this background depends on the solar cycle, and as the Sun is more 
quiet now, we expect a slightly larger astrophysical gamma-ray background at the present Fermi result than the older 
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EGRET measurement. 

Either way, comparing these background fluxes with the expected signal fluxes in Table HI we see that the expected 
signal is indeed, not only very small, but also extremely well hidden in the backgrounds. The proton-induced back- 
ground is expected to be quite well-defined to the solar rim, whereas the inverse-Compton up-scattered photons are 
expected to be more extended. In both cases, the Sun itself should appear as a shadow as it blocks the cosmic-rays 
from the other side. The WIMP signal on the other hand, should arise also from the Sun itself as is seen in Fig. \E\ 
However, even if we expect a different angular distribution and energy spectrum of the signal, it is going to be virtually 
impossible to disentangle the WIMP signal from the backgrounds, even for an overwhelmingly huge future gamma-ray 
detector looking towards the Sun. 



VI. DISCUSSION 



This work is under the assumption that the WIMPs in the solar halo are not disturbed during the capture process. 
In reality the Solar System also has planets, whose gravitational interaction can disturb the orbiting WIMPs. This 
effect, as caused by Jupiter, has been discussed in e.g. [l^. A massive planet, like Jupiter, can disturb orbiting 
WIMPs so that they end up on orbits no longer intersecting the Sun, potentially making them very long lived. On the 
other hand, WIMPs on orbits stretching so far out spend almost all their time far away from the Sun and hence do 
not contribute so much to the density close to the Sun. Furthermore, planetary interactions arc more likely to throw 
WIMPs out of the Solar System, reducing the WIMP density. WIMPs disturbed to bound orbits not intersecting 
the Sun end up on orbits sensitive to planetary interactions and will, in time, be disturbed again and typically lose 
their future contributions to the density close to the Sun. There is, however, a possibility that some fraction of these 
WIMPs on non-solar-crossing orbits turn out to be more stable, rendering a possible source for higher WIMP densities 
in the Solar System. Such a mechanism has been discussed in [s^, HJ , where it was argued that WIMPs which scatter 
in the outskirts of the Sun can be perturbed by the planets to orbits not crossing the Sun, in such a way that the 
new orbits become stable over very long timescales. This WIMP population would then build up over time and it 
was argued in [s^, HH to be a significant source of WIMPs in the Solar System. These calculations were analytical, 
but they have not been confirmed in later numerical simulations [s^l- In fact, in [32], it is argued that the life-times 
would be much shorter and that the density enhancement would be quite small. In principle, it would be interesting 
to add planetary interactions to our simulations, as these could in principle increase the density somewhat (even if 
not as much as originally argued in [sol. Isij). That is however outside the scope of this paper, as it would be a very 
extensive task to do. It is also unlikely that the effects would be large enough to be interesting. 

Another possibility to get higher WIMP densities is from a dark matter disk [33] , which could enhance the capture 
rate in the Sun by up to an order of magnitude (SJ) . This is predominantly because WIMPs from the dark disk pass 
at lower velocities and are hence easier for the Sun to capture as the energy loss required in the scatter is lower. 
However, lower initial velocities also imply that the just captured WIMPs will typically have lower energies, altering 
the distribution in Fig. [51 This implies that the captured WIMPs will survive less scatters before being fully trapped 
inside the Sun, which could reduce the density contribution per captured WIMP. Hence, the dark disk scenario need 
not increase the density around the Sun as much as it increases the capture rate. Even under the optimistic assumption 
that the density is effected as much as the capture rate, which would increase the gamma-ray signal by up to two 
orders of magnitude, the signal is still too low to be detectable. 

The calculated gamma-ray flux at Earth from WIMP annihilations around the Sun is extremely low, as seen in 
table m The calculated flux, at best, corresponds to around one gamma photon per (10 km)^ per century, which is 
hardly ever detectable. 

The low calculated signal in this work is in contradiction with the earlier result published by Strausz [9]. The 
analysis conducted here is, however, far more detailed. In it is assumed that the WIMPs move in one-dimensional 
orbits and that they all have some sort of average properties, such as the average energy loss in a scatter and the 
average scatter probability in a passage through the Sun. It is, however, very unlikely for the assumptions in Q to, 
by themselves, be responsible for such a strong disagreement. The passage where most of the disagreements appear 
is not explained in detail in and hence, it is impossible to tell what the source is of the difference in results. 

One should note also, that a calculation of the WIMP densities in the halo has also been performed by Fleysher 
[l^ . Fleysher flnds even larger densities than Strausz, but that calculation cannot be correct. The solar halo WIMP 
densities in the calculation by Fleysher are proportional to the scattering cross sections, which is unreasonable, as the 
scattering cross sections that enter in the probability of the flrst scatter and the lifetime of the WIMP orbits cancel 
(see section HV Al for a more thorough discussion about this). 

The conclusion of the work presented here, on the other hand, agrees with the conclusion argued by Hooper [ll|. 
That work is also not very detailed but more recent than Strausz' paper. Also Peter [2^ has performed simulations 
of WIMPs in the Solar System, from which the WIMP density can be extracted. As discussed in those results 
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are very similar to the results found here. Another test of the vahdity of our resuhs is performed in Appendix \X\ 

Finally, we can note that we set out this work with the hope of proving Strausz' optimistic calculations to be correct, 
but unfortunately, our detailed calculations instead show that the fluxes are indeed far too small to be detectable. 
On top of that, there are astrophysical backgrounds in which the signal is very well hidden. To reach this conclusion 
we have made optimistic assumptions about the WIMP dark matter candidate. However, we have assumed that the 
WIMPs are thermally produced in the early Universe, i.e. has an annihilation cross section of the order of 3 x lO^^e q^^s 
s~^. Of course, non-thermal productions or other enhancements (like a Sommerfeld enhancement at low annihilation 
velocities) would also be possible, but given the large enhancements needed to get an observable flux of gamma-rays 
from WIMP annihilations around the Sun, there are far better ways to search for those WIMPs in those cases, like 
gamma-rays from dwarf galaxies, annihilations in the Galactic halo, or effects on the cosmic microwave background 
radiation [3611. 



Appendix A: Are the results reasonable? 



The calculated signal is truly very weak. This appendix is a further check if the results are reasonable by also 
estimating the density using simpler arguments. Using what we have learned about the capture process and the shape 
of the WIMP density profile one can, assuming that these qualitative observations are valid, quite easily make an 
upper estimate of the magnitude of the WIMP density. This section studies the WIMP configuration M = 100 GeV, 
fSD = 10^"^ pb and (Tgi = 0. 

We can compare our total capture rate of WIMPs in the Sun with that estimated from previous work f20|. For 
the WIMP configuration chosen the WIMP capture rate is then Fc = 2.4 x 10^'^ s^^. Using this result, an upper 
estimate of the WIMP population available to build up the WIMP halo can be found by estimating an upper limit of 
the typical time scale for the WIMP capture process. 

From the Monte Carlo simulations one finds it to be extremely unlikely for a WIMP of this kind to need more than 
50 scatters for complete solar entrapment, which is easily verified to be reasonable by looking at the kinematics of 
the scatters. We assume that our hypothetical WIMP will lose the same amount of energy in each of its 50 scatters, 
which gives an upper estimate of the entrapment time because the assumption boosts the proportion of large orbits. 

The next assumption is that our WIMP has zero total energy before its first scatter. By this assumption the largest 
bound orbits are not taken into account, however, these orbits do not contribute much to the WIMP density close to 
the Sun, since a particle in this kind of orbit spends most of its time in the outer regions of the orbit and moves at 
maximum speed in the inner part. 

The reduced energy of the WIMP needs to be less than ^ —10^^ m^s^^ for the orbit to be completely trapped 
inside the Sun. This should here happen at scatter number 50; the WlMP's reduced energy after scatter number k 
is then £{k) = —2 x 10^ k m^s"^. The WlMP's energy after its first scatter then gives an orbit stretching out to 
?'max = 95i?0, which is more exact than the estimate above since the orbit needs a comparatively very low angular 
momentum to still cross the Sun, and Tmax can then be determined quite well. 

From these assumptions one gets an upper estimate of the total time the entrapment process takes a WIMP 

_ 1 ttGMq ^ 1 _ 1 ttGMq ^ 1 ^ 7.7 X 10^ 

P V2 fr[ (-f 4 X 1027/2 /t3/2 P ^' ^ ' 

where P is a lower estimate of the average scatter probability per solar passage for a WIMP in a bound orbit. The 
time required for a revolution in a given orbit is given by Eq. (j20p . It might be worth noting that this total time for 
the capture process also includes the time the WIMP spends inside the Sun. 

The total number of WIMPs orbiting the Sun at any given moment is then N — FcT^ot ^^'^ discussed above 
they maximally stretch out to 95P0. Assuming that the shape of the WIMP density curve in Fig. is correct, the 
WIMP line density is proportional to r^^ "*^. This gives the WIMP number density 

0.527V 2.0 X 1023 ^^_3 ^^2^ 



47r(950-52 - I)p0j52 p 
and hence, the gamma-ray flux at Earth, cf. ([55)1 

7.7 X 10-25 



*7 = — ^52 m-2s-i. (A3) 



What is left to estimate is the average scatter probability per solar passage. In the Monte Carlo simulations it was 
found for our WIMP scenario that the vast majority of the orbits has a scatter probability between 10""^ and 10~^ 
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per solar passage, orbits with scatter probabilities of less than 10~^ are found to be very rare. In this approximation 
we set P = 10~^. This corresponds to the scatter probability for a WIMP that passes through the Sun on a straight 
line with minimum distance to the solar centre of 40% of the solar radius, which makes a reasonable lower estimate 
since most orbits come closer than that to the solar centre. It is easily verified that this probability is reasonable for 
the column density in the Sun. 

For this WIMP configuration the upper estimated fiux of gamma-rays more energetic than 1 GeV is at Earth then 
finally 

$^ = 8 X 10"" m"2s"\ (A4) 

which is as expected a few orders of magnitude higher than the calculated flux from the Monte Carlo simulations of 
8.4 X 10~^^ m~^s~^, cf. table m The upper limit on the gamma-ray flux is too low to be detected. This reasoning is 
valid as long as a few long-lived WIMPs do not dominate the WIMP density. 
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